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ABSTRACT 


The use of the digital computer for the simultaneous generation of 
circuit responses and their corresponding sensitivity to parameter changes 
is considered for circuits containing diodes and transistors. Modeling of 
the diodes and transistors is based upon the Ebers-Moll equations, with 
the addition of voltage-dependent capacitors to account for switching 
time. The resulting general iterative equations are then used to calcu- 
late the effects of radiation on a p-n junction. The waveforms and 
circuit recovery time are studied as a function of circuit parameters. 

The use of sensitivity functions for predicting waveform variations is 


considered. 
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I. INTRODUCTION 


A measure of the sensitivity of circuit-response variables to param- 
eter variations is useful in circuit analysis and design. Many different 
definitions of sensitivity are available. H. Bode [1] defines a logarith- 
mic or normalized sensitivity function of a variable, x;(t.a,), with 


xX. onl a 


saat J ; , 
respect to a parameter, ass as D5. amin oe . The inverse of this 


definition is adopted by I. M. Horowitz [2]. Calculation of the sensi- 
tivity of linear-circuit responses to incremental changes in network 
parameters has been presented by J. Leeds [3]. Recently S. R. Parker [4] 
generalized the technique to include nonlinear circuits. Both Leeds and 
Parker define the sensitivity function as 51 : a . This work was 
preceded by that of P. Kokotovic [5] and R. Tomovic [6] who considered 
variations in the solution of a linear differential equation due to incre- 


mental changes in the coefficients. Tomovic defines the sensitivity 


Xe 3X. ; x 
coefficients as a = al and Kokotovic adopts the definition S ~ 
Fo :; ; 
Tat and discusses implementation of the sensitivity coefficient 

J 


using analog computation. The work by Leeds and Parker is oriented 
towards the digital computer. In this thesis digital calculation of 
Sensitivity functions of nonlinear circuits is considered in detail with 


specific application to the study of the effects of radiation on solid- 
Ne 
State junctions. The sensitivity function is defined as 5 "= is 
OX. j j 
7 oe ; nS 
Cias/a,y This definition is adopted since it enables sensitivity func- 
tions with respect to parameters having widely different nominal values to 


be conveniently compared on a percentage basis, and also avoids some 
difficulties encountered in the definitions adopted by Bode and Horowitz 


when the circuit-response variable is zero. It should be noted that the 


1] 


above definitions of the sensitivity function apply only for infinites- 
imally small changes in Xs and as The practicality of this definition 
when applied to specific examples with large parameter increments is 
considered in this thesis. 

In Chapter II it is shown that for linear circuits the sensitivity 
function may be derived using the compensation theorem [7]. A convenient 
set of iterative equations is developed using the concept of state vari- 
ables which allows the time response of the states and their corresponding 
sensitivities to be obtained simultaneously. A linear example illustrates 
the procedure. 

In order to obtain useful results from the sophisticated circuit anal- 
ysis and design programs for nonlinear circuits, acceptable nonlinear 
models are required. Chapter III discusses modeling of diodes and tran- 
Sistors for digital computer analysis based on the Ebers-Mol1 [8] equations. 
These models contain nonlinear capacitances and nonlinear resistances for 
the p-n junction. A state-variable formulation is developed which results 
ina set of iterative equations whose solution yields time responses for 
circuits containing linear elements, diodes, and transistors. The proce- 
dure is illustrated by considering the effect of a radiation current pulse 
injected across (1) a p-n junction diode, and (2) the base-emitter p-n 
junction of a transistor operating in the active region. These examples 
yield new results concerning the effect of radiation on p-n junctions and 
comparative data for several runs is presented. 

Sensitivity analysis of nonlinear circuits is discussed in Chapter IV. 
The state-variable formulation developed in Chapter III is extended, re- 
sulting in a set of iterative equations which allow time responses and 


their corresponding sensitivities to be obtained simultaneously. These 


12 


procedures are applied to the diode example of the previous chapter and 
the results interpreted in terms of the sensitivity of the diode response 


to percentage parameter variations. 


Is: 


Ii. SENSITIVITY ANALYSIS OF LINEAR CIRCUITS 


A. INTRODUCTION 

There are several general theorems which are useful in analysis and 
design of linear networks. One of these, the compensation theorem [7], 
has been known for many years and may be applied when it is desired to 
determine what effect an incremental impedance change in one branch of a 
network has upon the currents and voltages of other branches of the net- 
work. Recently the concept of the sensitivity of certain circuit-response 
variables to incremental changes in given network parameters has been 
presented by J. Leeds [3]. This work was preceded by that of P. Kokotovic 
[5] and R. Tomovic [6] who considered variations in the solution of a 
linear differential equation due to incremental changes in the coeffi-. 
cients. Tomovic and Kokotovic implemented the sensitivity coefficient 
by using analog computation. Leeds’ work is oriented towards the digital 
computer. 

In the next section sensitivity analysis of linear circuits is con- 
sidered. In the succeeding section the relationship between the older 
compensation theorem and the recent sensitivity concept is discussed. 

This correlation has not appeared in the literature heretofore. The final 
section contains an illustrative example to demonstrate the implementation 


of the results on a digital computer. 


B. DETERMINATION OF SENSITIVITY FUNCTIONS 


1. Auxiliary Network Approach 


Leeds [3] developed a relationship between the sensitivity func- 
ie ne 
tion, S_| = Beeb (where X 3 1S a response variable and a. 1S a ClPeunE 


a - 0a - 
J J 
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parameter), and node voltages and branch currents of an auxiliary network 
called the sensitivity model. The auxiliary network is obtained by re- 
ducing all independent sources to zero and adding a single source in the 
branch containing the element under consideration. The source depends 
upon the nature of the element. Fig. 2.1 summarizes the sources to be 
added for branches containing resistance, inductance, and capacitance 
respectively. The voltages and currents in the auxiliary network are the 
voltage and current sensitivities with respect to the parameter being 
considered. 

The advantage of Leeds’ approach is that it allows the sensitivity of 
a network to be calculated by evaluating the responses of a coupled auxil- 
lary network rather than actually performing the differentiation indicated 
by the definition of the sensitivity function. 

2. Augmented State Equations 

Sensitivity functions may also be obtained using the concepts of 

State variables. The network is first represented in state-variable form 


as discussed by T. Bashkow [9], and Kuh and Rohrer [10]. 
x = Ax + Bu (2.1) 


where 
“es yanenevector of cirewmt states 
A is a constant nxn matrix 


Heisei) VeClOr Gmmnonerng, TUunCcTOnS. 
Xs OX « 
The sensitivity functions, S.. = ——— 
a; a In a 
performing the indicated partial differentiation on (2.1). The resulting 


» may be obtained by 


partial derivatives can be considered as additional states which may be 
solved in conjunction with the original state equations. They may be 


included in the following augmented matrix; 


ES 


ORIGINAL CIRCUIT BRANCH 


CN eee 


AUXILIARY CIRCUIT BRANCH 


R , le 
dle 
+ dt 
O—VWYY\-—_(_ 0 
G 
or 
dt 


FIG. 2.1. CIRCUIT PARAMETER BRANCHES AND 
SENSITIVITY MODEL BRANCHES. 








. at | [ j 2 
ee ie a B u 0 0 
| | 
° 0 [a ‘ 
75) ! Sy | Uey | As) B sy 
| re i 
ic i ies ees E Xx 
| ’ : i: ; ol i220 
Xo e | i | 
=o ; 7 “sj | 1 de5 | Oe yi “ 
0 : — | 
; ’ Hu B 
et Ay) Asr BY! Ue | Ay sr 
or in partitioned form 
x Au Oualex B' 0 [ju fe x | 
anal awl a] eee — = = + ee a (2.3a) 
e { i 
sj a = y, BD HsJ| : "sj si za 


where j = 1,2,...,r , and r is the number of sensitivity parameters and 





aX) Xn ; 
a ne oe 
J at 
a) 
Qu Qu 
“sj “|oMay °° Tha, ee 
i 
| 
Beh 3411 
1 5 Wita. on > In a. 
Pee) ee J J (2.3d) 
SJ a Ina. ; ; 
J ° 
047 38 
0 Ina. d Ina 
<) I 
; bay Pam : 
bo In a- In a. 
ieee .| | 2] 
SJ 9 Inas. | mY (2.3e) 
J 3b 3D ay 
at a. 9 In a5 | 


0 
oe 7 | (2.338) 
0 -° | 
A 
B 
= (253ay 
By = B 


If u 1s independent of a» that is if a, is a circuit parameter 


only, (2.3a) takes the form 


| x< 
eee 
\! 

> 

>) 

| >< 

<—) 


(2.4) 





Ken Me 
a Sj 


Separating the set of first-order partial differential equations 


of (2.3a) for each sensitivity parameter yields 


Xp AX + BU eee a 


fOr 3) =: Fee eee 
Trapezoidal integration is commonly used to solve a set of first- 
order matrix differential equations because the procedure always converges 
and provides acceptable accuracy for most problems [11]. Utilizing trape- 


zoidal integration, (2.1) and (2.5) yield the following iterative solutions: 


x(nT) = g(T)x((n-1)T) + r(T)[u(nT)+u((n-1)T) J 
Xo(nT) = o(T)xXo.((n-1)T) + r(T) Lu (nT )tug ((n-1)T) J (2.6) 
+ Py CT Ex(nT)+xC(n-1)T) T+ r55(T)Lu(nT)+u((n-1)T) J 


TOW JH) sc els (2.7) 
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where 


(T) = [1 -5Al +54) 

MT) = (I-A & 
ryg(T) = [1 - pal! “sj 
ros(T) = [I - 5 Al” “sj * 


Equations (2.6) and (2.7) are a set of iterative equations whose 
solution yields the time response of the states and their corresponding 
sensitivities. Combining the results into a single equation results ina 


solution format similar to (2.3). 





/x(nT) oT)! 0 Hxt(n-1)T]. | ir(t)! 0 ||(u(mt) + uf(n-1)T]) 
— = eee Se i. S — 4 a 
ag (mT) 0 op (T) X, j[(n-1)T] Lo Tp (7) (ug s(nT) + us s[(n-1) 

of 0 | | (x(nt) + xf(n-1) TI) 
ry 4 (T) | Pp5(T) (X 5. nT) ¥ Xoi (n-1)T]) (2.8a) 
for = | Ze iG 
where _ 
o(T) | 
" 0 | 
by(T) 7 i 4(T) (28D ) 
0 x 
o(T) | 
r(T) 
0 
Ve nT) (2.8c) 
| 0 : 
L i) 


C. USE OF THE COMPENSATION THEOREM TO DERIVE SENSITIVITY FUNCTIONS 

The compensation theorem relates an incremental impedance change in 
one branch of a network with corresponding incremental current and voltage 
changes in the other branches [7]. In order to relate the compensation 
theorem directly with sensitivity analysis an analagous relation for incre- 
mental changes in the branch parameters R, L, and C is developed. 

1. Resistive Change 

If the resistance in a branch is changed by an amount AR, and a 

voltage source equal to I-(AR) is introduced (Fig. 2.2a) the network vari- 
ables remain unchanged. If another source is introduced (Fig. 2.2b) the 
branch voltage and current will change by an incremental amount AV and al 
respectively. By applying the superposition theorem (Fig. 2.2c) the 


incremental voltage and current can be calculated. 
AV Wome CR + AR et 1 (AR) (2.9) 


Dividing (2.9) by aR and letting AR+O results in (2.10). 


ov = ae (2.10) 


2. Inductive Change 


If the inductance in a branch is changed by an amount AL, and a 
voltage source equal to AL - (ty is introduced (Fig. 2.3a) the network 
variables remain unchanged. If another source is introduced (Fig. 2.3b) 
the branch voltage and current will change by an incremental amount AV and 
AI respectively. By applying the superposition theorem (Fig. 2.3c) the 
incremental current and voltage can be calculated. 

dal 


AV = (L+ AL) are + al => (22s 





Acti ve 


Network 


T+AI R+YAR I+AlI L+AL 


Active 


Network 















FIG. 2.2 RESISTIVE CHANGES. FIG. 2.3 INDUCTIVE CHANGES. 
I I+AI 
Grae dv 
| 
(a) (b) 
Ley 
Passive 
Mac 2 


Network 


(co) 
FIG. 2.4 CAPACITIVE CHANGES. 
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Dividing by AL and letting AL+0 results tn (2.12). 


ol. d] 
(<r) ae 2. 12) 


3. Capacitive Change 
If the capacitance in a branch is changed by an amount AC, and a 
current source equal to nce (SE) is introduced (Fig. 2.4a) the network 
variables remain unchanged. Introducing another current source (Fig. 2.4b) 
results in incremental changes in the branch voltage and current. Again 
by applying the superposition theorem (Fig. 2.4c) the incremental voltage 


and current can be calculated. 


A] =U en) ave Ac * (2.13) 


Dividing (2.13) by aC and letting aC+0 results in (2.14). 


(eee (2.14) 


Equations (2.10), (2.12), and (2.14) are identical to the auxil- 
jary network equations of Fig. 2.1. This illustrates the direct 


relationship between the compensation theorem and sensitivity analysis. 


D. LINEAR EXAMPLE 
Consider the peaking circuit [12] in Fig. 2.5a operating in the linear 
region. The incremental equivalent circuit is represented in Fig. 2.5b 


where > is the plate resistance and C is the coil capacitance plus 


p 
output capacitance, stray capacitance and wiring capacitance. The equiv- 
RR 
alent circuit is indicated in Fig. 2.5c where R = ees and v= 
ee eR 
R, Ro ] 2 
ere UL Vea 
Ry 4 Ro 1 


ae 


Ro L 
9) { 
V; Yo 
| | 
o— 
V 


(q) 





(c) 


FIG. 2.5 PEAKING CIRCUIT. 
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It is of interest to determine how sensitive the output is with 
respect to incremental changes in R, L and C. As representative values 
let L = 0.1 henry, C = 250ufarads, and R be adjustable according to the 
type of response desired. The response can be determined from the differ- 
ential equations of the circuit and is underdamped for K< 1 (R = 20), 
critically damped for K = 1 (R = 10), and overdamped for K > 1 (R = 5); 
where K = op /L/C 

The sensitivity functions may be obtained by constructing auxiliary 
networks as outlined in the preceding section, or by differentiating (2.1) 
with respect to In R, In L, and In C respectively. 

The state-variable representation is illustrated below. The solution 


takes the form of (2.8) where 























u =v = input voltage of 70 volts (2.15a) 
xy output voltage 
x= = (2.15b) 
Xo inductor current 
-1/RC -1/C lene 
A = .— B= (2.15c) 
lee 0 0 
Ue? 0 for j = 1.2% (2.15d) 
a a i 
aX) | sensitivity of output with respect 
dinR LOM 
Xo | = (2.15e) 
oe sensitivity of inductor current with 
5 In R eres pects EOmE | 
a ima sensitivity of output with respect 
9 InL-: EO 
Aso - - | a | (2 Jae 
2 sensitivity of inductor current with 
9 Ink | respect Temi 
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sensitivity of output with respect 











a alin C ton C 
X53 = = | (2.159) 
aK) sensitivity of inductor current 
a In Ci] with respect to C 
r 7 7 7 - 
1/RC 0 i O 0) ake & 1/C 
Ag) * | > Agy =| leis Hae =i 
fom Gl ey oe on 0 
(2.1.50) 
-1/RC 0 | TARE) 
B = B = B = Zar) 
$] : S2 7 $3 | 
= 0 oe 


The computer program used to implement the solution is included in 
Appendix B. This program is general in that it is easily modified to 
handle linear circuits containing more state variables and sensitivity 
parameters. 

The time response of the output ie the circuit and the corresponding 
sensitivities are presented in Figs. 2.6 through 2.11 for underdamped, 
critically damped and overdamped responses. These results have been 
obtained using trapezoidal integration as well as a fourth-order Runge- 


Kutta integration technique. 
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FIG. 2.6. LINEAR EXAMPLE --UNDERDAMPED RESPONSE. 
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FIG. 2.7 SENSITIVITY FUNCTIONS FOR UNDERDAMPED 
RESPONSE. 
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FIG. 2.8. LINEAR EXAMPLE-- CRITICALLY DAMPED 
RESPONSE. 
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FIG. 2.9. SENSITIVITY FUNCTIONS FOR CRITICALLY 
DAMPED RESPONSE. 
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FIG. 210. LINEAR EXAMPLE --OVERDAMPED RESPONSE. 
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IIIl. MODELING OF NONLINEAR CIRCUITS 


A. INTRODUCTION 

The digital computer has played an increasingly important role in 
circuit analysis and design in the past ten years. This has been partly 
due to the emergence of sophisticated circuit-analysis and design programs 
[13]. In order for these programs to be useful, acceptable nonlinear 
models which relate model parameters to physical processes are required. 
These models must lend themselves to numerical analysis and provide accept- 
able quantitative accuracy. 

In the next section, a general mathematical format which is convenient 
for sensitivity analysis of nonlinear circuits is developed. In the suc- 
ceeding section modeling of diodes and transistors is outlined. This jis 
followed by examples which consider calculation of the response of a 


p-n junction to a pulse of radiation current. 


B. GENERAL NONLINEAR STATE-VARIABLE FORMAT 


The state equation for a nonlinear time-invariant system 1s of the form 


x(t) = a(x(t),u(t)) (Sai 


where a is a vector of nonlinear functions of the state variables and 
forcing functions. 

As is shown in the section which follows, (3.1) takes the following 
general form for circuits containing diodes and transistors (represented 
by the Ebers-Moll [8,14] equations). This representation includes non- 


linear resistances and nonlinear capacitances. 


x(t) = A(x(t))ex(t) + e(x(t)) + B(x(t))-u(t) (3.2) 
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where A(x(t)) is a nxn matrix and B(x(t)) is an nxm matrix, both with 
elements which may be nonlinear functions of the states, and c(x(t)) is 
a general nonlinear vector. 


Equation (3.2) may be written as 


x(t) = #(x(t)) + B(x(t))- u(t) (3.3) 


where 


|—h 
—, 
—, 


bh) = Bet -meet) + alxGt)) 


(x 


Using trapezoidal integration for (3.3) results in 


(n) = x(n-1) = 4D [e(x(n)) + £(x(n-1))] 


| >< 


+ SF [B(x(n))-u(m) + B(x(m-1))-u(n-1)] (3.4) 


Equation (3.4) is a nonlinear matrix equation implicit in x(n) which, 
as a result, is generally difficult to solve. Simplification results if 
f(x(n)) and B(x(n))-u(n) , where u(n) is considered constant, are expanded 


in a Taylor series about the point x(n-1). In matrix notation 


f(x(n)) = 5, (ax'o)* #(x(n-1))/k:! (3.5a) 
B(x(n))-u(n) = z (ax! a) k [B(x(n-1)).u(n) J/k! (3.5b) 
k= 
where 
3 3 i 
i Lox aX> x 


7, 10, X, WX are W vecters 
and 


me 4 XU - x(n-1) 


Bo 


A fimsit-order approximation for (3e5) resus man | oooe 
Fy) = FUT yp AX 
B(x(n))-u(n) = B(x(n-1))-u(n) + (ax 


Substituting (3.6) into (3.4) yields 


(i) Bess 1-11) An [ 2#(x(n-1)) + (ax! 


— 


V) f(x(n-1)) J 


|x 


+ a [B(x(n-1)) > (u(n) + u(n-1)) + (ax'v) {B(x(n-1))-u(n)}] 


(3 ge) 
Equation (3.7) may be simplified by using the matrix identity 
(ax'v) f# = (vf)! ax (3.8) 
which is proven in Appendix A. 
Substituting (3.8) into (3.7) and solving for x(n) results in 
eT ea ree -Attvel)! + {vfB-u(n)] 37! 
[f + 5B ° (u(n) + u(n-1))] (3.9) 


where f, B, vf! and (B-u(n))! are evaluated at the point x(n-1). 
Equation (3.9) can be solved directly for x(n) from the known values of 


x(n-1), u(n), and u(n-1). 


C. DIODE AND TRANSISTOR MODELING 
1. Modeling for a p-n Junction Diode 
Diode modeling has received much attention during the past twenty 
years [15,16]. Recently Steele [17] prepared a report concerning semi- 
conductor modeling specifically for computer analysis. The applicable 


results are summarized in this section. 
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a. Current-Voltage Relationship 


The current-voltage relationship for an ideal diode is 


; _ av 
io I, (e - 1) (3.10) 
where I, = reverse saturation current 
a = a parameter dependent upon the temperature and the 


type of semiconductor material. 


Generally the reverse saturation current is not entirely 
independent of voltage. This is due to a leakage current and can be 
accounted for by the addition of a resistance, Rois shunting the diode. 
Also, at large revérse valtages an effect known as Zener breakdown occurs.. 
This effect is omitted in the following analysis. 

b. Diode Resistance 

The static resistance Rtv) of a diode is defined as the ratio 
V/i 4 Since Rylv) varies widely with v and igo it is generally not a 
useful parameter. An incremental, or dynamic, resistance rv) is defined 


as rylv) : ow = . Although the dynamic resistance varies 
d 


Ioae™ 





widely with v it is a useful parameter for small-signal operation. 
c. Capacitance Effect 

Diode capacitance is due to two separate effects; space charge 
and charge storage. 

(1) Space Charge. A diffusion potential, Vp? is established 
in the junction region when p and n materials are placed in contact. This: 
1s caused by an increase in donor atoms on the n side and acceptor atoms 
on the p side of the junction, resulting in a capacitive effect. As a 
reverse voltage is applied more donor and acceptor atoms are uncovered 
causing the junction width to expand, and hence the capacitance to decrease. 


This voltage-dependent capacitance has the following form when the diode 


oD 


is reverse biased. For radiation analysis, acceptable results are ob- 


tained when Cy is assumed to be constant. 


| o 
Cy = Ga (35am 
where n = 1/2 for abrupt junctions 
n = 1/3 for linear-graded junctions. 


(2) Charge Storage. When the diode junction is forward biased 
the potential barrier is lowered. This results in a high density of minor- 
ity carriers injected into the p and n regions, which can be viewed as 
stored carriers. A change in voltage results in a change in density (num- 
ber of stored carriers) and hence a capacitive effect. 

If the voltage on the junction is abruptly reversed the 
stored carriers reach equilibrium after some storage time Ty The associ- 


ated current may be defined as 


pee (3.12) 


Substituting (3.10) into (3.12) yields 


a - av | dv 
i = Thi ae = C.(v) Te (Soler 
where C.(v) - Tae is the charge-storage capacitance. The charge- 


storage capacitance may be alternately written in terms of the dynamic 


resistance as C.(v) = Ty/r gv) 


d. Diode Model 
Combining the effects outlined above results in the diode 


model in Fig. 3.1. The circuit equation becomes 


mn dv os V 
C4(v) ag oe (3.14) 


where C yy) = Ch 7 C(v}. 


<= 


Fig. 3.1. Diode model. 


oy 


Cp 


C(v) 


Equation (3.14) is in the form of (3.3) where 


cae) (3.15a) 


I 
| 
) 
is 
<= 
©) 
Cy 
<= 
7O 
—Y 
a 


Bxt ts) ew (3.15b) 


The diode example (III. D.) does not include parameter temper- 
ature dependence. The range of values studied includes a spread in 
parameter values as may be due to manufacturing tolerances, environmental 
conditions, and aging. The Zener breakdown effect and voltage dependence 


of C, are neglected. Typical parameter values are used. 


D 
2. Modeling for a p-n Junction Transistor 


Many models for a p-n junction transistor have been proposed in 
the last fifteen years, the most common being the Ebers-Moll model [8,14], 
the lumped model [16], and the charge-control model [18,19]. It has been 
shown that each model yields the same result in the solution of a large- 
Signal transient problem [20]. The Ebers-Moll model is used in this study 
and is summarized in this section. 

a. Current-Voltage Relationship 

The current-voltage equations for the Ebers-Moll model are 


given in (3.16) and (3.17). 


a7! I 
I = po (eC 1) - ~—EP— (eME - 1) (3.16) 
Ob Orn O17 On 
ont EO Ov C0 Ov 
I. = a (e aoe 1) = lara ae (e C = Wy (301s 
OT ON TON 
where 
a = constant having a value between 1 and 2 


= common-base current in the normal mode 
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= common-base current gain in the inverse mode 


I = saturation current of collector junction with zero 
emitter current 


I = saturation current of emitter junction with zero 
collector current 


b. Transistor Model 
Figs. 3.2a and 3.3a are representations of (3.16) and (3.17) 
for a n-p-n and p-n-p transistor respectively. Figs. 3.2b and 3.3b repre- 
sent the corresponding equivalent circuits if each diode is modeled as 
outlined in the preceding section. 


The reverse saturation currents are as follows: 


a,l 
rac 
I = (3.18a) 
Onleg 
ES 1 - a.a 
Pew 


If the positive senses of the currents are as indicated in 


Page. oe and 3.3, leo and Leo are positive quantities for n-p-n transistors 


dj 
e e e s a dc 
and negative quantities for p-n-p transistors and Coo (v) = The 7" 
— di : ® ° ° 
and Coe (v) = The — . The is approximated by the transistor rise 


time and T.,. iS approximated by the sum of the transistor storage time and 


DC 
the transistor fall time. The rise, fall and storage times are complex 
functions of the base current, collector current, Cys Oly» and the transistor 
current-gain cut-off frequency, and are defined in Fig. 3.4 [17]. 

If Cav) = Che + Con (v) and Cue (Vv) = Coe + Coe (Vv) the 


equations for Fiq. 3.2b become 





= 0 (3.19a) 


o 
es 
oO 
~~ 
< 
— 
er 
ctr 
-- 
—! 
Bs 
O 
+ 
A 


a9 





(a ) 





de | de 
(oa hata 
way Cy 
DOC DE 
ic Cy. (v) Cop (v) ce 
A L. E 
lp 
Za (_ 
a : oy) 
ldc B 
R suc 
(b) 
vc > ec 


Fig. 3.2. Model and equivalent circuit for the n-p-n 
transistor. 
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Co lv) I. 








iG. ag (b ) oe i 


= 


Fig. 3.5. Model and equivalent circuit: for the p-n-p 
transistor. 
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lp = base current 


Time 


(a ) 


l.= collector 
current 





rise storage fall 
Time 


Fig. 3.4. Transistor switching times. 
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=. 


V V 





E E c 
Cue (V) ct tse + - Shige. = = 0 (3.19b) 
SHE 
which are in the form of (3.3) where 
zs 
ON'de de “C \ 
Cal) Cac!) Roucoage™) 


(1 'de ae fe 
Cael) Cael) RcHECde lV) | 


B a [0] (3. 20b) 


Similar equations can be written for Fig. (3.3b) 
The transistor example (III. E.) does not include the effects 


of parameter temperature dependence and Cy» C T,- awd Tp are assumed 


be” 
to be constant. Typical parameter values are used. 


DE C 


D. DIODE EXAMPLE 

A diode problem of considerable interest is now solved in order to 
illustrate the modeling and the mathematical formulation discussed in this 
chapter. The problem considers the effect of a radiation current pulse 
injected across a p-n junction as indicated in Fig. 3.5, where the radiation 
current pulse is defined by 


di. (it) 


ai +e CU UE! (3.21) 


where y(t) is the unit pulse given in Fig. 3.6. 
Replacing the diode in Fig. 3.5 by the model described above yields 


Fig. 3.7. Typical values are 


I. = 1.0 pA Cy = 120005 Ty = ].0us 
J = a z 

a = 1/.026 V Rou = 100.0 MQ Tp = Q.2us 

A= 10.0 mA 
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R=lOk 7\t) 


F= lOv lpp | 
| time (us) 
Fig. 3.5. Diode circuit. Fig. 3.6. Radiation current 
pulse. 
R 
: Cy Colv) 
WL = Q) 
a V ly 'pp 


Fig. 3.7. Equivalent circuit. 
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Analytical Solution for Radiation Current 





Equation (3.21) can be solved using Laplace transforms. 


| | ~/ 
A:( — - ——~~— ) 0<t< 10 
I (s) = ; oP i (3.22) 
pp 1 Sy, 
394 Ar( TT t > 10 
The resulting time response is given in (3.23). 
e/ T, OF a t= 1077 
( A-(1 -e R) 
ie = (Sco) 
~/] 
omnes 7 t> 107 


Piecewise~Linear Diode Model 


If the diode current-voltage relationship is assumed to be piecewise 


linear (Fig. 3.8) ig =JOMfomay—0vandiy = igRy for v > 0, where Rie Rylv) 


ry(v). The storage capacitance becomes a constant, C. . "D/R , for 


v > 0 and an open circuit for v < 0. Fig. 3.7 can be redrawn as in Fig. 


3.9 where Ry is the parallel combination of R, Rou and Ra and C = Co fier 


V 


> 0, and Re is the parallel combination of R and Row and C = Cy for Vv <0. 


The rise and fall time can be determined using Laplace transforms. 


For 0 <t < 10°/ 
ion ee (3. 24a) 
s + T7RC 
a E/R 
ICS) = I p's? = ee (3.24b) 
aca lees eg eran 


S oa I/Tp S 


The voltage takes the form 


7 


v(t) = C, + Coe VIR + Cae S/R One t <a wereees) 
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Fig. 3.8. Piecewise linear model. 





Fig. 3.9. Diode circuit representation. 
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For t > 1077 


: (alo Ba : 
i .394 A.(— 7. é (3.26) 


The resulting voltage response is given by 


v(t) = C, + cae /TR + es t a te? (3.27) 
The differential equation becomes 
ee = it a er (3.28) 
dt pp . 


where G = 1/Ry. 


Equation (3.28) is solved using trapezoidal integration. Several 
values of Ry are used and the output voltage, Vi = is indicated in 
Fig. 3.10. The computer program is presented in Appendix B. The solution 
in Fig. 3.10 is of little practical value because of the wide variation of 
the voltage response with assumed values of Rye 


3. Taylor Series Approximation 


Equation (3.28) may be rewritten as 
Coyidv = wt. dt - Givi-vwdt - E/R dé (3.29) 


where C(v) = C, + C.(v) 


G(v) == + =— + 
R Roy Rg Cv) 


[c(v(n)) + C(vin-1))] J = So Ciyain) + iy(n-1) - I 
(3.30) 
- SF [e(vin))-v(n) + G(v(n-1))-v(n-1)] 
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OUTPUT VOLTAGE (volts) 





TIME (ps) 


FIG. 3.10. NONLINEAR EXAMPLE -- PIECEWISE LINEAR 
MODEL WITH Rg ADJUSTED FROM 
O1kKQ TO 4.0 kQ. 
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Expanding C(v(n)) and G(v(n))v(n) in a Taylor series and including 


only the first-order terms yields 


oC 


Geum = —C(v(n-1)}) + 2 eae wy = C(v(n-1)) (3.31a) 
G(v(n))-v(n) = G(v(n-1))-v(n-1) + G(v(n-1))-Av + 2 re li 
vin- 
= G(v(n-1))-v(n) (3.31b) 
where it is assumed Bh encdy QO and 5 hina) OF. 


Substituting (3.31) into (3.30) results in 


c(v(n-1))-av = SO Ci (n) + aaa (ner) - 2 
i (3.32) 
- 5 HONS) tv) --v(n-1))] 
Solving for v(n) yields 
an) = 2C(v(n-1)) - AT-G(v(n-1) ) me 
2C(v(n-1)) + AT-G(v(n-1)) 
iin) See ne) = Se 
— ‘OOD R (5.33) 


2C(v(n-1)) + AT ° G(v(n-1)) 


It is of interest to note that the same result is obtained if 
(3.28) is integrated over the interval (n-1)T to nT as a linear, constant- 
coefficient equation with C(v) and G(v) having the values C(v(n-1)) and 
G(v(n-1)) respectively. 

This method is implemented with an adjustable integration step 
size. The step size is adjusted to insure that C(v) and G(v) do not vary 
by more than a fixed percentage from (n-1)T to nT. The resulting output 


voltage, Was is indicated in Fig. 3.11, where the maximum that C(v) and 
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OUTPUT VOLTAGE (volts) 





FIG. 3.11. NONLINEAR EXAMPLE --DIODE RESISTANCE AND 
CAPACITANCE ARE HELD CONSTANT BETWEEN 


CALCULATION POINTS. PERCENTAGE VARIATION 
INDICATES THE MAXIMUM ALLOWED VARIATION 


IN C(v) AND G(lv) BETWEEN CALCULATION POINTS. 
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G(v) are allowed to vary is 10%, 1% and 0.1%. The computer program is 
presented in Appendix B. This method becomes more accurate as the percent- 
age variation allowed for C(v) and G(v) is reduced. 
4. Diode Solution using General Nonlinear State-Variable Format 
The circuit in Fig. 3.7 contains a single state variable, v(t). 


The state equation for the circuit can be represented in the form of (3.3) 


where 
- G-v(t) - I, (gov (t) - ® 
f(v(t)) = | aa (3.34a) 
Cy sr Tplae 
] 
B(v(t)) = a (3.34b) 
Ch “ TT .ae 
u(t) = Toe - E/R (3.34c) 
] ] 
wammesG = — + 
Ese 


For the single state-variable problem considered in this example 


(3.9) reduces to 


aT [f +B * (u(n) + u(n-1))] 


muy = — _ AT pat, a(B uta) —) 
2 “OV OV 
where the partial derivatives may be obtained from (3.34) and f, B, a, 


and Ss are evaluated at v(n-1). 


The computer program for this method appears in Appendix B and the 
results are indicated in Fig. 3.12. The program adjusts AT automatically 
so that Av never exceeds 0.01 volts. These results have been confirmed 
by direct integration of the nonlinear state equation using a fourth-order 


Runge-Kutta integration scheme. 


2 | 


OUTPUT VOLTAGE (volts) 





TIME (ys) 


FIG. 3.12. NONLINEAR EXAMPLE -—- SOLUTION USING 
TAYLOR SERIES METHOD. 
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5. Presentation of Data 

The results of several runs are presented in Figs. 3.13 through 
3.17. The recovery time, defined as the time for the voltage to recover 
to 90% of its original value, is plotted as the ordinate. The abscissa 
represents the parameter value. 

In Fig. 3.13 the amplitude of the radiation pulse, A, is varied 
from 0.265 mA to 20.0 mA. When A < 0.265mA the radiation current pulse 
does not cause the voltage to drop below the 90% value. 

inkigs. 3.14 Whr@egh 3.17 Cys R, Tp and Tp are adjusted with 
A = 5.0 mA and A = 10.0 mA respectively while the other circuit parameters 


are held constant at their nominal values. 


E. TRANSISTOR EXAMPLE 

The effect of a radiation current pulse injected across the base- 
emitter p-n junction of a transistor operating in the active region (Fig. 
3.18) is considered in this section. The radiation current pulse is 
defined by (3.21) where y(t) is the unit pulse given in Fig. 3.19. 

1. Transistor Solution using General Nonlinear State-Variable Format 


Modeling the transistor in Fig. 3.18 by the equivalent circuit 


presented in Fig. 3.2 results in Fig. 3.20 where 


avec 
wae) ~ Anees@e (3.36a) 
- ee ae (3.36b) 
eee = - 1d (3.36c) 
“de = Ipe (2 © - wp (3. 36d) 
Typical parameter values used are 
Oy = | lee = 0.01nA Toe = |].0nsec 
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Recovery time (ys ) 





OF 
Cy (pF) 
Fig. 3.14. Diode recovery time as a function of Cp . 
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Fig. 3.15. Diode recovery time as a function of R. 
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Fig. 3.16. Diode recovery time as a function of T,. 
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Fig. 3.17. Diode recovery time as a function of Tp . 
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Collector 





Emitter 


Fig. 3.18. Transistor circuit. 


r (t) 


1.5 
time (ys) 


Fig. 3.19. Radiation current pulse. 
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Emitter 


Fig. 3.20. Equivalent transistor circuit. 


60 


a, = 0.99 Ice = 0.1pA ieeO as 


Tp = 0.4 nsec CoE = 1.0pF a = 38.46 V 
Roc = 100 MO Che = Te0pF 
Roue = 100 MQ A = 0.1mA 


The circuit equations are obtained by summing the currents at the 


collector and emitter nodes yielding 








ae ee 
Cuelv) At tee + ue Op de 7 cae 0 (3.37a) 
SHE 
dv Ve Ve 
C,(v) =—+ +i, -a,i, +i, +— = 0 3.37b) 
dc t Rowe dic N de R Rp 
where 
] 
ip = ne Wee +. Ve = Ve) (3.37c) 


The biasing resistors, R, and Ros are calculated for a desired 
biasing voltage vane base=collector voltage Ve; and base-emitter voltage 
dv dv 
ees ; E ; 
Ve under steady-state conditions, by setting iS saci and ‘yp 


equal to zero in (3.37). The program used is presented in Appendix B, 
and the resulting values are R, = 0.35 M2 and R. = 4.5 kQ for Nee = 10.0 V, 
0.6 V and Ve = -3.4 \V. 

Substituting (3.37c) into (3.37a) and (3.37b) and solving for 


os 


x = [ve Ve v yields a result identical in form to (3.3) with 


G 
Gr: | 
Cn) C, (v)- 
f(x(t)) = de" an anaes 
ete 
| Catv) "Re Cue (v) 


6] 





av ae OVE 7 
(= lic (e° Cuaap) ; Ipc (e 1) 
( (v) c (v) 
* = oS (3.38a) 
QV a OVE 2 
( ote Seay I.p(e 1) 
Sancen 7 
Cac : Cac(Y) | 
Wy 1 
Caan V Cte aR 
B(x(t)) =| °° os 
0 ] (3. SSb» 
Cac (Vy) R 
1 o 
pp! 
u(t) = (3.38c) 
| Ue 
] ] ] ] ] 
where G =— 1+ =— and G See a eee 
ie Ro Rou ine Re Ro Rouc 


The iterative solution equation is of the form of (3.9) where 


ve! and v(Bu(n))! are obtained by performing the indicated gradient opera- 


EAOMmONsies oe 
The computer program for this problem appears in Appendix B, and 
the results are indicated in Fig. 3.21. The program adjusts AT automatically 


so that the norm, |e] , never exceeds 0.01 volts. 


ia /(ovg)” + (Ave)* (3.39) 

These results have been confirmed by direct integration of the 
nonlinear state equations using a fourth-order Runge-Kutta integration 
scheme. 
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OUTPUT VOLTAGE (volts) 


> 


N Gl 
a 


0.0 
0.0 


FIG 3.2. 


TIME (us ) 


TRANSISTOR RESPONSE TO RADIATION 
CURRENT PULSE. 
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2. Presentation of Data 

The results of several runs are presented in Figs. 3.22 through 
3.25. The recovery time is plotted as the ordinate and the abscissa 
represents the parameter value. 

In Fig. 3.22 the amplitude of the radiation pulse, A, is varied 
from zero to 3.0mA. When A < 2.0uA the radiation current pulse does not 
cause the voltage to drop below the 90% value. 

InebidsSees.23 threvah 3.25 The: The: Coes Cor and Tp are varied 
respectively while the other circuit parameters are held constant at their 


nominal values. 
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Recovery time (us) 








loci mio’ sec) 
Fig. 3.23. Transistor recovery time as a function of Tog and Toc. 
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Recovery time ( ws ) 


Recovery time (us) 





Fig. 3.24. Transistor recovery time as a function of Cpe and Coc 


sO} 5. 10. 15. 20. 25. 


Tr (nsec 


Fig. 3.25. Transistor recovery time as a function of Tp - 


IV. SENSITIVITY ANALYSIS OF NONLINEAR CIRCUITS 


A. INTRODUCTION 

Sensitivity analysis of nonlinear circuits and systems is a relatively 
new field. Recently S. R. Parker [4] has considered the sensitivity analysis 
of general nonlinear circuits utilizing auxiliary coupled networks. The 
sensitivity analysis of nonlinear circuits utilizing the concept of state 
Space is covered in the next section, followed by a sensitivity analysis of 


the nonlinear diode circuit discussed in Chapter III. 


B. SENSITIVITY-AUGMENTED STATE EQUATION 
The state equation format adopted in Chapter III (repeated below) is 


convenient when considering circuit sensitivity. 


x(t) = A(x(t))-x(t) + e(x(t)) + B(x(t))-u(t) (4.1) 


xe dX. 
As in the linear case, the sensitivity function, Se = Wet ; 
J J 
may be obtained by performing the indicated partial differentiation on 


(4.1), yielding the sensitivity-augmented state equation. 





x A'o {lx c PB, OMNlu 0 ee 
aes air 2 1+ j2— |e Sete ee + |l eee 
Xsi} [2 1 “oy lSsi} [Ssa} 1° Poff 4sa} Asa, Psi] ] 4 
(4.23 
where j = 1,2, ....¥r and X oq Ans By» Aga Beye Ug. are defined in Chapter 
II and 7 
—_ et 2 eee oan (4.3) 
= S$} d In a. d In a. a Tn a. 


Separating the set of first-order partial differential equations for 


each sensitivity parameter in (4.2) yields 
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Xoa = A Xs. pace As; X sR Us, + BS u (4.4) 


ner yp =*8.2, ..9,F 


Equation (4.4) may be rewritten as 


Xsj e Fo g(XsX55) + B(x)u.. + By (XsXea)u (4.5) 


where j = 1,2, ... ,r and 


f .(%3X0.) = A x _. ew Tee y 
=S) > s 5) oreo — 


Applying trapezoidal integration to (4.5) yields 


Ke g(M) = Xes(n-1) = FE [Fe (x(n) x55(m)) + Fo5(x(n-1),x,5(n-1))] 


+ SF (B(x(n))ugs(n) + B(x(n-1) )ugs(m-1)] + FY (B,5(x(n) x65(n)) ulm) 


fa Be. (X(N 1) x 


53 .(n-1)) u(n-1)] (4.6) 


*sj 


The functions £4 (x(n) x<.(n)) and B (x(n) x¢5(n))u(n), where u(n) 


is considered a constant, can be expanded in a Taylor series about the 


point (x(n-1),x..(n-1)) and approximated by the first two terms, as in 
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Chapter III, yielding 





f 4 (x(n) .x¢5(n)) = fc. (x(n- 1) X54(n-1)) 
+ (ax" 9, )F 5 (X x(n-1) ,x X¢,(n-1)) + (axe 964) fg 5 {Xm-T) X55 (n-1)) (4.7a) 
Be j(x(n) .x<-(n))u(n) = Bo .(x(n-1),x.-(n-1) )u(n) 
+ (axT9,) [By <(x(m-1) x 5(n-1))u(n) J + (0x2 505 )1B G5 (x(n-1) 9x, (n-1)) ul) ] 
(4.7b) 
where 
= i 
7 = Se a 
, 
a 3 
xsj : ox aXe | 
SJ SJ 


B(x(n))-u. s(n), where U. 5 i) 1s considered constant, can be expanded 
in a Taylor series about = point x(n-1) and approximated by the first 


two terms yielding: 


B(x(n))-ue s(n) = B(x(n-1))-u.s(n) + (dx'¥, B(x(n-1))-u, s(n) J 


=a) =o) a 
(457i 
Substituting (4.7) into (4.6) yields 
we Al il, 
x, 3 (n) 7 X 3 (n 1) 2 L2f,. 5 * (dx ile i (ax ,v, 
Al F 
+ => [B (u, 5 (n) oan 5 (n Tax" ve )(Bu, .(n))] 
+ Al re Muth) Fe utn-1)) + (ax voy(B_.u(n)) + (ax) fy.) (Bue 
Ss) = — Xx! sSje —sj xsj’\ sj- 
(4.8) 
where ite ee and Be ju(n) are evaluated at the point (x(n-1)), 
Recall that 
(Ax! Tee <ul ee ee ae AX, (4.9) 


sd Kes) 55 


Substituting (4.9) into (4.8) and solving for X. 5 (n) yields 


eae: ai alY * aT [I - A ik eng Z (Vg 0B, ;u(n) ]!)"} ie 


(ax! 9 ) fi aa (u, 5 (n) + U. je 1)) + (Ax! ve ) (Bu, j(n)J 
+B. -(u(n) + u(n-1)) + (ax'v,)(B, ,u(n) D3 J (4.10) 


Equations (3.9) and (4.10) can be solved directly for x(n) and x. (n) 


from the known values of x(n-1), u(n), u(n-1), u.j(n), and u u.(n-1). The 
solution yields the time response of the states and their corresponding 
sensitivities. 


/0 


C. NONLINEAR EXAMPLE 

In order to illustrate the technique discussed above, sensitivity 
analysis of the diode circuit in Chapter III is considered. The state 
equation for the diode circuit of Fig. 3.5 contains a single state var- 


jable. Equation (4.10) reduces to 





lias 
v_.(n) = v_.(n-1) + AT“ Lf. + 7 >t Ay + Be (u, .(n) 7 u.;(n-1)) 








SJ SJ OV 
a(Bru. n) 3(B. ,u(n)) 
7 el Av + B57 (u(n) + u(n-1)) + - -— Av} ] 
one 3(B_.u(n)) 
mT Sj Sj -| 
(l- = { +—=1____ }}, J (4.11) 
2 OV 5 OV | 
ae 3(B-u_.(n)) 6(B_.u(n)) 
= Sc 2J | 
where ie a > B, Beas By and =a we evaluated 


at the point v(n-1). 


The state equation may be written in the form of (4.1) where 


A(v(t)) = - —+__ (4.12a) 
Cy + Tpigae 
[ce - aah) 

c(v(t)) = - 2———_ (4.12b) 
Ch + Tyl ae 


and B(v(t)) and u(t) are given in (3.34b) and (3.34c) respectively. 
The equations for the sensitivity functions take the form of (4.5). 
The sensitivity parameters considered are the diode parameters, Ty and 


C., the radiation-pulse parameters, A and Tp» and the circuit parameters, 


D? 
Rand E. The terms in (4.5) are presented in table 4.1 for each of the 


sensitivity parameters where 


Tak 


7] 
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The expressions for ead rans are obtained by performing 





the indicated partial differentiation on (3.21) and = ; Dy : 

9(B_.u(n)) of 3(B_ .u(n)) 

ee ; and _ are obtained by performing the 
OV OV 6 | a j 


indicated partial differentiation on the equations listed in table 4.1. 
In addition to the sensitivity functions the voltage response, 

= TN is shown for (1) nominal parameter values, and (2) a 10% change 

in the sensitivity parameter value. The response is presented in Figs. 

4.2, 4.4, 4.6, 4.8, 4.10, and 4.12. The computer program for generating 

the sensitivity functions is given in Appendix B, and the resulting 

sensitivity functions are presented in Figs. 4.1, 4.3, 4.5, 4.7, 4.9, and 


4.11. The program adjusts AT automatically so that je, | never exceeds 


0.05 volts and |Av| never exceeds 0.01 volts, where 


i 


= /(ay)* + (av, )° Eat j=) ee... AIG (4.13) 


These results have been confirmed by direct integration of the non- 
linear state equations using a fourth-order Runge-Kutta integration scheme. 
It should be noted that the method outlined in this chapter is approxi- 


mately twice as fast as the fourth-order Runge-Kutta integration scheme. 


D. INTERPRETATION AND APPLICATION OF THE SENSITIVITY FUNCTION 


The definition of the sensitivity function adopted in this thesis is 


x. OX. OX. 
i 1 1 Sa itas 
a; = om a, 5a, 7a, - As noted, this definition applies only for 


infinitesimally small changes in X. and as. The sensitivity function may 
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FIG. 4.1. SENSITIVITY FUNCTION WITH RESPECT TO Tp. 
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FIG. 4.2. RESPONSE OF DIODE CIRCUIT TO RADIATION 
PULSE. 
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4.5. SENSITIVITY FUNCTION WITH RESPECT TO E. 
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FIG. 4.4. RESPONSE OF DIODE CIRCUIT TO RADIATION 
PULSE. 
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FIG. 4.5. SENSITIVITY FUNCTION WITH RESPECT TO R. 
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FIG. 4.6. RESPONSE OF DIODE CIRCUIT TO RADIATION 
PULSE. 
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FIG. 4.7 SENSITIVITY FUNCTION WITH RESPECT 
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FIG. 4.8. RESPONSE OF DIODE CIRCUIT TO RADIATION 
PULSE. 
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FIG. 49. SENSITIVITY FUNCTION WITH RESPECT TO Cp. 
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FIG. 4.10. RESPONSE OF DIODE CIRCUIT TO RADIATION 
PULSE. 


83 


SENSITIVITY FUNCTION (volts) 
Oo Oo S 


= 10) Pt fo ff 


TIME (ys) 


FIG. 4.11. SENSITIVITY FUNCTION WITH RESPECT TO Tp. 
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FIG. 4.12. RESPONSE OF DIODE CIRCUIT TO RADIATION 
PULSE. 
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be used to approximate changes in a response due to fractional changes 


in a parameter as follows 


xe Aa 


- 7 j 
AX. = S.°(—) (4.14) 
7 a a 


In general this approximation becomes more accurate as the percentage 
change in the parameter value approaches zero. Table 4.2 represents a 
comparison of results obtained for the underdamped response of the linear 
example in Chapter II. The output voltages, Vy 2 V7 2 +25 9 Ve » are 
obtained using resistance values of Ro? la) Ro? MERRANGS Ro? where Ke 
is the nominal value. Approximations for these voltages are calculated 


using (4.14); specifically, 


V. = Vv eS ee) (4.15) 


~ 


where vee is the approximate response after Ry has been changed by ARG» 

and Vo is the nominal response. Table 4.3 presents a similar comparison 
for the nonlinear diode example. Values are tabulated for time intervals 
when the sensitivity function does not vary rapidly. The results in 

tables 4.2 and 4.3 indicate that the approximation for the change in output 
voltage improves as the amplitude of the sensitivity function approaches a 
constant value. In this region correct values are predicted for variations 
as large as 50%. It is noted that (4.14) does not yield useful results 

for parameter variations as small as 1% when the sensitivity function dis- 
plays sharp peaks. In this region the sensitivity function may still be 
used to determine changes in circuit responses by assuming that the asso- 
ciated parameter variations cause a time shift rather than a voltage 
change. Fig. 4.13a represents an idealized switching-circuit response 
where the effect of an incremental parameter variation is upon the switch- 
ing time. Fig. 4.13b represents the corresponding sensitivity function. 
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go 


Bau 8.9] 
2.0 9.37 
1.0 9.76 
Sensitivity 
Function 
OV 
OTe Vo 
(volts) 
-10.0 9.02 
-10.0 9.24 
-10.0 9.83 
Table 4.3. 


8.59 8.61 8.27 8.31 7.90 7.4] 
9.14 9.17 8.91 8.97 8.19 8.37 
9.65 9.66 9.52 9.56 9.06 "omen 
ne AE ae 
E I E c cn * 
Ly Va Vo Ps Mis V4 
10.03 10.02 11.03 11.02 14.03 14.02 
10.24 10.24 11.24 11.24 14.24 14.24 
10.83 10.83 11.83 11.83 14.83 14uee 


Voltage comparison for nonlinear diode 


example using (1) actual parameter changes, (2) approximation 
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FIG. 4.13. (a) IDEALIZED CIRCUIT RESPONSE 
AND (b) ASSOCIATED SENSITIVITY 
FUNCTION. 
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The delay time, Ty due to the parameter variation is given as 
Ty = t - t, (4.16) 
Applying (4.14) to Fig. 4.13 it follows that 


at 
ea eos = h har ) (45g 


The area, A, under the sensitivity function over the time interval 
ie 


ma t, 1S 


X. 
si) dt = hTy (4.18) 


Solving for h in (4.17) and substituting into (4.18) yields an 


equation for the delay time uae 


(Aa s/a,) ‘ 
le = Saas” f > 5 dt (4.19) 
Ss t, J 


Equation (4.19) may be applied to estimate the delay time correspond- 
ing to a pulse-type sensitivity function by calculating the area of the 
pulse and substituting into (4.19). This type of calculation is performed 
for the waveforms of Figures 4.1, 4.3, 4.5, and 4.7 for 10%, 20% and 50% 


parameter variations. The results are compared with the actual time 
delays in table 4.4. The area, A= f ; si 
that the amplitude spikes are eee Times ty and ty are chosen at 
those points where 7 = .15 h . The results indicate time delays can be 


dt , is calculated assuming 


calculated within 10% for parameter variation as large as 50%. It is 
noted that calculation of the area presents the greatest limitation in the 


accuracy of the delay time. 
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V. CONCLUSIONS 


The technique of generating sensitivity functions for linear circuits 
using the concept of a sensitivity model is shown to be an extension of 
the compensation theorem. A set of iterative equations is developed to 
yield simultaneously the time responses of the circuit states and of their 
corresponding sensitivities. An example illustrates the procedure. 

Diode and transistor modeling for discrete computation based upon 
the Ebers-Moll equations is considered. The models adopted contain both 
nonlinear capacitors and nonlinear resistors for the junctions. A state- 
variable formulation is developed which results in a set of nonlinear 
iterative equations for circuits containing linear elements, diodes, and 
transistors. The use of these equations is illustrated by considering the 
effects of a radiation pulse on a p-n junction, and the resulting recovery 
times of the junction are tabulated. 

The foregoing nonlinear iterative equations are extended to generate 
simultaneously the sensitivity functions as well as the circuit solution. 
The preceding example, which considers the effects of a radiation pulse on 
a p-n junction, is then analyzed from the point of view of a sensitivity 
analysis. This analysis indicates that proper interpretation of the sen- 
sitivity function results in accurate estimated for changes in circuit 
responses due to parameter variations as large as 50%. When the amplitude 
of the sensitivity function is constant the change in the circuit response 
due to a parameter variation is reflected as a change in response ampli- 
tude. However, in regions where the sensitivity function displays impulse- 
like amplitude peaks, the change in the circuit response due to a parameter 


variation is reflected as a time delay or shift in the response characteristic. 
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APPENDIX A 


PROOF OF (ax'y) ¢ = (vel)? ax 
The matrix equations ax! vf may be written as 
T : 3 = 3 
(Ax Vv) i: = [ Ax, Eau T AX» aX + AX Some 
] Z n 
| ei af af, | 
ON ey OnaK, Xn OK 
| af, fF of, | 
Axe .— + ie, —— AX, a=) 
] OX, 2 9X5 N ox, 
af of of | 
n n n | 
Dike. Se AX, = + AX 
] OX] Z OX » n An 
Consider this matrix equation 
~ oof, of 
T | clin eae 
ao =e [ft fd enor ea 
OX, | : at We 
| ee 
| es 
| as 
ax 
n i 
i 3%, «OX y 


oo 


ae 


(A.1) 


(A.2) 


‘Taking the transpose of equation(A. 


of. 


TT a 





AX) 


AX 


of 


Z 
OX 
n 


AX» ae 


AX» wae 


of 


n 


AX Sere 
2 Xo 


2)and multiplying by Ax yields 


AX 


(A.3) 


Equating equation(A.3) to equation(A.1) results in equation(A.4), which 


completes the proof. 


iti 
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APPENDIX B 
COMPUTER PROGRAMS 
5 UTION AND SENSITIVITY ANALYSIS PROGRAM FOR LINEAR 
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